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The phase behavior of hehcal packings of thermoresponsive microspheres inside glass capillaries is 
studied as a function of volume fraction. Stable packings with long-range orient ational order appear 
to evolve abruptly to disordered states as particle volume fraction is reduced, consistent with recent 
hard sphere simulations. We quantify this transition using correlations and susceptibilities of the 
orientational order parameter '06- The emergence of coexisting metastable packings, as well as 
coexisting ordered and disordered states, is also observed. These findings support the notion of 
phase transition- like behavior in quasi- ID systems. 



The phenomenology of ordered phases and phase trans- 
formations in systems with low dimensionality is surpris- 
ingly rich. While dense three-dimensional (3D) thermal 
packings can exhibit long-range order, for example, this 
trait is absent in purely one dimensional systems [1]. 
Complexities arise, however, when considering 3D sys- 
tems confined primarily to one dimension (ID). Investi- 
gation of order and phase behavior in quasi- ID thermal 
systems, therefore, holds potential to elucidate a vari- 
ety of novel physical processes that have analogies with 
polymer folding [2] , formation of supermolecular fibers in 
gels [3] , and emergence of helical nanofilaments in achiral 
bent-core liquid crystals [4]. 

Packings of soft colloidal spheres in cylinders provide 
a useful model experimental system to quantitatively in- 
vestigate order and phase transformations in quasi- ID. 
At high densities, spheres in cylindrical confinement are 
predicted to form helical crystalline structures [5,6]. Ev- 
idence for such packings have been found in foams [7,8], 
biological microstructures [5], colloids in micro-channels 
[9] and fullerenes in nanotubes [10]. However, research 
on these systems has been limited to static snapshots and 
athermal media. Recent simulations suggest that transi- 
tions between different helical ordered states [11,12] and 
between quasi-lD ordered and disordered states [13-15] 
should exist in thermal systems, but such transitions have 
not been investigated experimentally. 

In this Letter, we explore ordered and disordered struc- 
tures in a quasi- ID thermal system of soft particles with 
adjustable volume fraction. In particular, we create heli- 
cal packings of thermoresponsive colloid particles in glass 
microcapillaries, we show theoretically that phases with 
long-range orientational order can exist in quasi- ID, we 
demonstrate experimentally that such phases with long- 
range orientational order exist at volume fractions below 
maximal packing, and we study volume-fraction induced 
melting of these orientationally ordered phases into liquid 
phases. The orientational order parameters and suscepti- 
bilities that characterize these phases and this crossover 
are measured and analyzed. Coexisting regions of or- 
dered and disordered states and coexisting ordered do- 



mains with different pitch and chirality are observed at 
these crossover points. Such coexistence effects suggest 
the presence of abrupt or discontinuous volume-fraction 
driven transitions in quasi- ID structures. Interestingly, 
these orientationally ordered phases in quasi- ID share 
physical features with orientationally ordered phases ob- 
served [16,17] and predicted [18] in 2D. 

The experiments employed aqueous suspensions of 
rhodamine-labeled poly-N-isopropylacrylamide (NIPA) 
microgel spheres (polydispersity < 3%) with diameters 
which decrease linearly and reversibly with increasing 
temperature [19]. The unique thermoresponsive charac- 
teristics of NIPA microgel particles provide a means to 
explore the phase behavior [17,20-21] of soft spheres in 
quasi- ID as a function of volume fraction. Borosilicate 
glass tubes (McMaster Carr) were heated and stretched 
to form microcapillaries with inner diameters compara- 
ble to NIPA microsphere diameter. An aqueous suspen- 
sion of NIPA microspheres was then drawn into the cap- 
illaries; subsequently, the capillary ends were sealed with 
epoxy and attached to a glass microscope slide. Sam- 
ples were annealed at 28° C to permit uneven packings to 
re-arrange at low volume fraction, thereby creating sta- 
ble high- volume fraction packings when returned to lower 
temperatures. 

The samples were imaged under a 100 x oil-immersion 
objective (N.A. = 1.4) using spinning-disk confocal mi- 
croscopy (QLC-100, Visitech, International). Resulting 
images depict 75/im long segments of densely packed re- 
gions which span at least several hundred microns. An 
objective heater attached to the microscope (Bioptechs) 
permitted control of the sample temperature to within 
±0.1° C. Standard particle tracking routines [22] were em- 
ployed to identify particle positions from three dimen- 
sional image stacks and 2D cross-sections. 

At high densities, we observe crystalline helical pack- 
ings with varying pitch and chirality dependent upon 
particle- and tube-diameter. The observation of large or- 
dered domains is consistent with the tendency for these 
nearly monodisperse particles to form uniform crystals in 
2D [17] and 3D [20]. Particles in hard-sphere helical pack- 
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FIG. 1: (color online) (a) Geometrical characteristics of ob- 
served structures, (2,2,4) (blue □), (1,3,4) (purple o), (2,3,5) 
(green A), (1,4,5) (cyan y), (0,5,5) (red O), (3,3,6) (dark 
green ★). Filled symbols indicate structures with observable 
Brownian motion. Dashed vertical lines indicate theoretical 
D/dsep values for predicted structures. Inset: Cartoon of 
axial cross- section of a (336) structure with dsep > ^eff' 
Confocal flourescence images of helical NIPA packings at high 
volume fraction. Scale bar = 10/im. 



ings predicted for cylinders [5] have six nearest neighbors 
whose relative order along the tube axis corresponds to 
a characteristic set of three integers (m, n, m + n). This 
notation for distinct helical crystalline structures is com- 
monly used in phyllotaxy and is similar to the vector 
used to describe carbon nanotube chirality. We verify 
such crystalline ordering from analysis of 3D confocal 
image stacks. All varieties of predicted helical packings 
in the given range of tube-diameter-to-particle-diameter 
ratios were observed in 15 samples (Figure 1), with the 
exception of structure (0,4,4). Additionally, the ratio 
D/dsep^ where D is twice the average radial distance 
from the particle centers to the central axis of the tube 
and dsep is the average nearest neighbor particle sep- 
aration, fell within experimental error of the predicted 
maximally packed hard- sphere values. 

Ordered structures were found to exist over a range 
of volume fractions below maximal packing. When 
the effective particle diameter d^jj^ [23] was such that 
dsep/d^jj < 0.95, the particles did not appear to move 
(i.e., motions greater than 0.2 jiia were not observed dur- 
ing the 10-second scan). In such cases, we consider the 
particles to be packed closer than their effective diam- 
eters. For dsep/dgff > 0.95, particles fluctuate signif- 



icantly about their equilibrium positions; thus, thermal 
helical crystalline structures exist at volume fractions be- 
low close packing. 

The volume fractions of such samples were then low- 
ered further to determine if, when, and how the pack- 
ings disorder to isotropic states. Two uniformly packed 
samples were chosen and are presented here for careful 
analysis: a (2,2,4) packing of microspheres with diame- 
ter 1.71/im at 22°C and a (0,6,6) packing of microspheres 
with diameter 1.23/im at 22°C. The sample temperature 
was increased in steps of 0.2-0.7°C. At each tempera- 
ture step, after allotting ample time for the sample to 
reach thermal equilibrium (at least 5 minutes), videos of 
two-dimensional cross-sections of the packings were taken 
at 15 - 30 frames per second for approximately 5 min- 
utes. Though these two-dimensional videos lose some of 
the structural information available in three dimensional 
scans, they provide data at much higher speeds and yield 
better axial position tracking of particles in view. 

A local orient ational order parameter, ipQj = 
^^'^'^ e^'^^^^ /Nnnj quantifies helical order in these sys- 
tems. Here, Ojk is the angle between the axis of the 
tube and the bond between particles j and /c, and Nnn 
is the number of nearest neighbors of particle j. Though 
this order parameter is typically used for two-dimensional 
planar systems, it is acceptable to use in the analysis of 
two-dimensional slices of helical packings. Helical pack- 
ings are effectively two-dimensional triangular lattices 
wrapped into cylinders, and the observed cross-sections 
of these particular packings exhibit only slight variation 
from the ideal two-dimensional hexagonal lattice. 

We examined the spatial extent of orientational order 
along the tube by calculating the orientational spatial 
correlation function, gQ{z = \zj — Zk\) = {tl^Qjil^Gk) ^ where 
Zj is the axial position of particle j. As depicted in Figure 
2, the resulting correlation functions decrease quickly at 
low volume fractions, as expected in a disordered state. 
However, at high volume fractions, these functions ex- 
hibit long-range order within the experiment's field of 
view. These experimental observations are consistent 
with an expectation of long-range orientational order. 
Though long-range translational order is impossible in 
one dimension at finite temperature [1], long-range ori- 
entational order is possible, just as in the much-storied 
theory of two-dimensional melting [24]. 

One can show theoretically that long-range orienta- 
tional order persists even in this quasi-one-dimensional 
system by evaluating the orientational correlation func- 
tion ^6(1") using the isotropic elasticity free energy [24,25] 
in an "unwrapping" of the particles on the cylinder sur- 
face onto a two-dimensional infinitely long strip [26]. We 
assume that particles are packed densely enough so that 
fluctuations along the radial direction of the cylinder may 
be neglected. As r ^ oo, ^6 approaches a constant. 
Thus, finite correlations exist at infinite distance, a hall- 
mark of a phase with long-range order. Orientational 
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FIG. 2: (color online) Local maxima of orientational spatial correlation function g6{z) for samples with (2,2,4) packing (a) and 
(0,6,6) packing (b). Dashed lines: Full correlation function at ^1 't>jnax(224) ~ ^'^^ ^'^^ ^/ ^max(066) ~ ^'^^ 

0.36 (b), where (jymax is the maximum volume fraction for a given packing of hard spheres with the observed particle spacing. 
Empty symbols indicate samples with short-range order. Oscillations arise from the periodic distribution of particle positions 
along the tube axis. Full correlation functions were calculated for all volume fractions, but omitted for the sake of clarity. 



order arises here through the crystalline axes defined by 
unwrapping the cylinder, rather than through an explicit 
additional mode, as in such systems studied in [27]. 

To our knowledge, the existence of long-range orienta- 
tional order has not been characterized in previous stud- 
ies of packing in cylindrical systems [11-15]. At low vol- 
ume fractions, one expects long-range orientational corre- 
lations to disappear as observed in the experiment. How- 
ever suggestive this combination of theory and experi- 
ment may be, we emphasize that further work is required 
to elucidate whether a generalization of KTHNY theory 
[24] is appropriate for this system. 

To quantify the crossover from long-range to short- 
range order in the experimental system, an average ori- 
entational order parameter is defined for each frame, 
where V^e = ;^ XI j V^e j • The height of the first peak 
of the one-dimensional axial structure factor S{kmax)-, 
where S{k) = [N{N - 1)]-^ |^ ^ e^^l^'-^-l |, is used as 
a translational order parameter. Here, z denotes the ax- 
ial position of each particle, N is the number of particles 
in the field of view at a given time, and kmax is chosen 
iteratively for each volume fraction. 

In Fig. 3(a,b), it is evident that both the average trans- 
lational and orientational order parameters cross from 
an ordered state at high volume fraction to a disordered 
state at low volume fraction, though the change in 
is significantly sharper than the change in Sikmax 

). The 

orientational crossover complements a recent simulation 
which finds a similiar crossover in hard sphere packings 
with decreasing density [14]. 

To characterize the fluctuations in orientational or- 
der, we calculate the orientational susceptibility xe = 
(iV^eP) — (iV^el)^, where () represents the time average. 
Statistical effects of the finite size of the system are ac- 
counted for by calculating the susceptibility of different 




FIG. 3: (color online) Top: Mean values of translational 
(red T) and orientational (blue ■) order parameters for (a) 
(2,2,4) and (b) (0,6,6) samples. Bottom: Orientational sus- 
ceptibilities for (c) (2,2,4) and (d) (0,6,6) systems, for sample 
size = 75 /im (o) and = oo (•). Dashed vertical lines at 

^l't'max(224) ^ ^'^^ ^"^^ I max (06 6) = 0-38. 



sub-segments of length L in the system and extrapolating 
to the limit L ^ oo^ similiar to the calculations in [17]. 
Plots of xe in Fig. 3(c,d) clearly demonstrate a peak in 
the orientational susceptibility. The location of this peak 
coincides with the onset of both orientational and trans- 
lational order in the system. Translational susceptibili- 
ties were also calculated, but did not exhibit clear peaks 
or trends with respect to the order parameters or volume 
fraction. We do not expect any transition-like behavior 
from translational susceptibilities due to arguments from 
[1] against long-range translational order, which apply in 
quasi- ID. 
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FIG. 4: (color online) Snapshots of particle tracks from sys- 
tems exhibiting emergence of coexisting domains, (a) (0,6,6) 
structure with emerging disordered domains. Color repre- 
sents local variance of the phase of V^e, (^ij) — (^6j)^, in 1.2 
/xm segments, (b) (2,3,5) structure with emerging domains 
of (3,2,5) and (0,4,4) structures. Color represents phase of 
V^Gj, ^6j, which characterizes packing orientation. Structures 
corresponding to each Oqj given on right. 



The existence of a diverging peak in the susceptibil- 
ity of an order parameter typically indicates a phase 
transition [28]. However, it is difficult from the given 
data points to determine if this is truly a diverging 
peak, and whether it would indicate a first-order (asym- 
metrically diverging) or second-order (symmetrically di- 
verging) phase transition. Upon closer examination of 
the (0,6,6) sample, we observe coexistence of ordered 
and disordered domains for 61 (i)^^rr,fnaa\ — 0.47-0.40 



max(066) 

(see Figure 4(a)). We also observe the appearance of 
a small domain with dubious order in the (2,2,4) sample 
at 4^1 4^fji(ix(224) ~ appearance of these coex- 

isting domains is consistent with the spatial correlation 
functions exhibiting neither long-range nor short-range 
behavior at intermediate volume fractions in Figure 2 

('^/<^maa:r066;=0-40' 'i'/'i'max(224) = O'^^)' '^^^ P^^^" 
ence of such solid-liquid coexistence has also been seen 
in recent simulations of hard spheres in cylinders [15]. 

In other samples, domains with different helical or- 
der often appear as volume fraction decreased (Figure 
4b). The appearance of these domains was difficult to 
quantify, since domains would grow, shrink and/or dis- 
appear with decreasing volume fraction. The structures 
observed in coexisting states were those with most similar 
predicted linear densities and D / dsep values, consistent 
with recent hard sphere simulations [11,12]. This coexis- 
tence of ordered structures should not be confused with 
dislocation-mediated structural transitions theoretically 
studied in [6] and observed in [8] in athermal helical crys- 
tals, especially because of the difference in the observed 
sequence of coexisting structures. 

In summary, we created ordered helical packings of 
thermoresponsive colloids and observed the presence of 



long-range order resilient to thermal fluctuations. Sharp 
crossovers from orientationally ordered to disordered 
phases with decreasing volume fraction were observed. 
In addition, we flnd basic evidence for abrupt volume- 
fraction driven structure-to-structure transitions. These 
flndings raise and elucidate fundamental questions on the 
subject of melting in ID. 
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HELICAL PACKINGS AND PHASE TRANSFORMATIONS OF SOFT SPHERES IN CYLINDERS: 

SUPPLEMENTAL MATERIAL 



NIPA MICROSPHERE DIAMETER CHARACTERIZATION 



Dilute suspensions of the same NIPA particles used in the tube packings were placed between glass coverslips such 
that the gap between coverslips was approximately the diameter of the NIPA particles, creating a quasi-2d monolayer. 
Videos of these particles diffusing in two dimensions were taken using a 100 x oil-immersion objective (N.A. = 1.4) at 
five temperatures from 24 to 28 °C. Particle centers were tracked using standard particle tracking routines [1]. The 
two dimensional pair correlation function g{r) of the particle locations was then calculated from the particle tracks. 
At each temperature, the approximate diameter of the particles taken to be the first value where g{r) = 1/e, since, 
in the first approximation, g{r) ^-u(r)/kBT ^ effictive diameter of particle is often taken as the value of r for 

which U{r) = IksT. 

Since previous studies have observed a linear relationship between diameter and temperature for these particles in 
this temperature range [2], we then take a linear fit of these data points to find a functional relationship between 
effective particle diameter and temperature. For the larger species used in this experiment, we find the relationship 
deff = 2.52/im-0.037(/im/°C) x T, and for the smaller species, we find 4// = 2.41/im-0.054(/im/°C) x T. 

For videos of two dimensional cross sections of particle packings in cylinders, the nearest neighbor particle spacing 
dsep was calculated from axial particle spacings and theoretical helical packing values. The volume fraction ratios of 
the tube packings were then calculated as (l)/(j)max = {def f / dsepT • 




r (microns) r (microns) 

FIG. 5: (a) Two-dimensional correlation functions for (a) smaller and (b) larger NIPA particles at different temperatures. 
Horizontal dashed lines at g{r) — 1/e. 
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FIG. 6: Approximate diameters for the smaller (blue •) and larger (red k) NIPA particles at different temperatures, with linear 
fits. 
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IDENTIFICATION OF STRUCTURES WITH PREDICTED HELICAL PACKINGS 



We describe the predicted helical packings with a commonly used three- index notation, (m, n, m + n) [3]. If we 
consider any single particle in a such an ideal close-packed system of hard spheres, we notice that each particle six 
nearest neighbors. We can thus think about such a packing as a two-dimensional triangular lattice wrapped around a 
cylinder. The indices then indicate the unit vectors in the unwrapped triangular lattice connecting any point to itself 
in the helical structure. Alternatively, one notices that the three indices indicate the relative distances of a particle's 
nearest neighbors along the axis of the cylinder. For example, if any given particle's nearest neighbors are the second, 
third or fifth closest particles in the axial direction, it would exist in a (2,3,5) packing. 

After identifying particle centers in 3D using common particle tracking routines [1], the nearest neighbors of each 
particle in an image are identified as those closer than the far end of the first maximum in the 3D pair correlation 
function g{r) (see Figure 3b). After an individual particle in the packing is selected, all other particle are given integer 
values based on the order of their axial distance from the selected particle in each direction. The integer values of 
the nearest neighbors are then recorded. This process is repeated for every particle in the structure. This creates 
histograms of the relative axial order of the neighboring particles. The peaks in these histograms then identify the 
integers (m, n, m + n) describing the ideal packing (see Figure 3c). 
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FIG. 7: Steps involved in categorizing structures and quantifying structural parameters of a sample (2,3,5) structure, (a) Three- 
dimensional reconstruction of microsphere packing from tracking of confocal image stack, (b) 3D pair correlation function from 
particle positions. First peak in g(r) highlighted in red, with Lorentzian peak fit in blue dashed line. Vertical dotted line 
indicates nearest neighbor cutoff, (c) Axial and azimuthal positions of tracked particles in a short axial section. Selected 
particle position highlighted in green, with nearest neighbors highlighted in blue. Relative axial distance from the selected 
particle listed to side. Note that the nearest neighbors have relative axial order of 2, 3 and 5. (d) Histogram of relative axial 
order for 1st, 2nd and 3rd closest nearest neighbors in the axial direction. Peaks at values of 2, 3 and 5 indicate that this is a 
(2,3,5) packing. 



ESTIMATION OF STRUCTURAL PARAMETERS FROM 3D TRACKS 



The average inter-particle spacing dgep in a 3D track is determined from the 3D pair correlation function. Since 
each structure has on the order of 100 particles, the resolution of said correlation function is fairly poor. Thus, we 
determine dsep by fitting a Lorentz peak, /(r) = /o + A{w/{4:{r — ro)^ + w'^)) to the first peak in g{r) (see Figure 3b). 
The center of this function, ro, gives us the value for dsep- The value for structure diameter D was taken as twice the 
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average value of the radial positions of the particles in the structure. 



ELASTIC APPROACH TO LONG RANGE ORIENTATIONAL ORDER IN A CYLINDRICAL SYSTEM 

OF COLLOIDS 



In one- and two-dimensional systems of particles, there can be no long range positional order, however, there can 
be long range orientational order in two dimensions [4]. This fact is key to the KTHNY theory [5-7] of continuous 
melting in two dimensions. 

We point out here that the same arguments used by Nelson and Halperin in [6] applied to a quasi- one-dimensional 
system, show that there is long range orientational order in such systems at low temperatures. The system we consider 
is a two-dimensional box which is infinite in the x direction and periodic with length L in the y direction. Note that 
such a system system applies to the problem of spherical particles in a cylinder if we assume that at high densities, 
the outer layer of spheres in the cylinder behave as if they were disks lying in a two-dimensional strip with periodic 
boundary conditions - that is, we "unwrap" the spheres onto a plane and neglect the radial (which become out-of- 
plane) fluctuations. This maps ijjQ deflned in the quasi-one-dimensional system approximately onto the ijjQ that was 
measured on the surface of the cylinder. 

We'll calculate a few quantities in the low-temperature ("solid" phase) by using an isotropic Lame elasticity free 
energy. The following expression is appropriate for the elasticity of a triangular lattice with the given geometry: 



F = - dx I dy [Xu% + 2ijiUijUij] 



(1) 



Here Uij = ^ {diUj + djUi) is the displacement from the (zero-temperature, hexagonal) ordered state. This free 
energy assumes that there exist nonzero elastic moduli A and //, a natural assumption we must make. 

To calculate the correlation function of orientational order, we first define ipe = e^**^, where measures the bond 
angles between nearest neighbor particles (the factor of 6 arises from the triangular symmetry of the lattice). Fur- 
thermore, we will use the fact that the angle can be shown to be |V x i/. 

The fact that the theory is quadratic allows us to take two shortcuts. First, we can evaluate {2pQ{r)ipQ{0)) (where 
r = (x^y)) as the exponential of an average. Second, this average can be evaluated simply since the inverse of the 



elasticity dynamical matrix is known [6,8] to be D-^ = \ [^jn — ) 



}■ 



It may be worthwhile to point out that we can guess the result from "counting powers of g". We can see that 
translational order is destroyed in this system by estimating the scaling of the (Fourier-space) correlation function of 
u at low q. To do this, we integrate DJ^ in one dimension. But DJ^ scales as so the result / q~'^dq scales like 
1/q. This diverges at low g, meaning that long- wavelength fluctuations drive positional correlations to zero. 

However, since 6 is related to the derivative of the (Fourier-space) correlation function of will involve an extra 
factor of q for each = \/ xu. We thus expect the two-point correlation function to scale as / dq{q^){q~'^) ^ g, which 
is flnite at low g, indicating long range order at inflnity. 

We now proceed with a more detailed calculation of the real-space orientational order correlation function: 



(V^6(r)V^6(0)) = exp[-36((^^(0)) - (^(r)^(O)))] 



= exp 



exp 



exp 



exp 



9kBT 
' 27rL ^ 



ij^mn / ^ I dq^qiq-m (1 ^ 



{qxX+2'Kny / L) 



n=-N 
N .A 



I dq^ (l - e^(9-^+27rn2//L)^ 

AT ^ — A 



N , 

J2 (2A-e'2''"^/^- 



'.Ax ^—iAx 



QkBT f ^ sin(27r(iV+l/2)t//L)2sin(A:r) 

2tt fiL \ sm{7ry/L) x 



(2) 



We now let the cutoff A be 27r{N -\-l/2)/L and use the sine function (sine t = (sint)/t): 
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('06(r)'06(O)) = exp 
= exp 



^ksT [ ^2 _ p2^y sin(A7/) sin(Ax) 



TT^/i \ L sm{7ry /L)Ay Ax 

9kBTA^ / sinc(Ax)sinc(A7/) Y 
TT^/i \ smc{7ry/L) J 



(3) 



For large this approaches the constant exp hence this quasi-one dimensional system has long 

range orient at ional order. This result explains our experimentally observed long range correlations in the orientational 
order parameter at high densities. 

We emphasize that (3) should not be fit to experimental or simulation data in the present form for the following 
reasons. First, the elastic field theory describes the sphere systems at large wavelengths. This approach is fine for 
the purposes of searching for long-range orientational order, but is inappropriate for the extraction of quantitative 
correlation functions. 

Another issue is that the precise functional form of (3) arises from the cutoff scheme chosen. Ideally, one would 
apply a physical theory which describes the higher wavelength modes and how the system responds to them (in the 
case of spheres, part of this scheme could be a density functional- like theory), which would be more suitable for 
capturing the behavior of the correlation function over short distances than the hard cutoff that we used. 

Thus the primary experimentally accessible result for comparison to theory is the fact that the correlation function 
remains finite at infinity. 
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